climdex.pcic (version 0.2-1)

climdex.su, climdex.fd, climdex.id, climdex.tr: Frost Days, Summer Days, Icing Days, Tropical Nights

Description

These functions compute the climdex indices FD, SU, ID, and TR.

Usage

climdex.fd(ci)
climdex.su(ci)
climdex.id(ci)
climdex.tr(ci)

Arguments

ci

Object of type climdexInput.

Value

• A vector containing the number of days meeting the criteria of the index for each year.

Details

These functions take a climdexInput object as input and compute the climdex indices FD, SU, ID, and TR. These are:

• FD:

{Number of frost days; annual count of days days where daily minimum temperature drops below 0 degrees Celsius.} SU:{Number of summer days; annual count of days where daily maximum temperature exceeds 25 degrees Celsius.} ID:{Number of icing days; annual count of days where daily maximum temperature is below 0 degrees Celsius.} TR:{Number of tropical nights; annual count of days where daily minimum temperature exceeds 20 degrees Celsius.}

References

http://cccma.seos.uvic.ca/ETCCDMI/list_27_indices.shtml

See Also

climdexInput.raw, climdexInput.csv.

Examples

## Create a climdexInput object from some data already loaded in and
## ready to go.

## Parse the dates into PCICt.
tmax.dates <- as.PCICt(do.call(paste, ec.1018935.tmax[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
tmin.dates <- as.PCICt(do.call(paste, ec.1018935.tmin[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
prec.dates <- as.PCICt(do.call(paste, ec.1018935.prec[,c("year",
"jday")]), format="%Y %j", cal="gregorian")

## Load the data in.
ci <- climdexInput.raw(ec.1018935.tmax$MAX_TEMP,
ec.1018935.tmin$MIN_TEMP, ec.1018935.prec$ONE_DAY_PRECIPITATION,
tmax.dates, tmin.dates, prec.dates, base.range=c(1971, 2000))

## Create an annual timeseries of the number of frost days.
fd <- climdex.fd(ci)

## Create an annual timeseries of the number of summer days.
su <- climdex.su(ci)

## Create an annual timeseries of the number of icing days.
id <- climdex.id(ci)

## Create an annual timeseries of the number of tropical nights.
tr <- climdex.tr(ci)

## Print these out for testing purposes.
fd
su
id
tr